Cosmological expansion and local physics 
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The interplay between cosmological expansion and local attraction in a gravitationally bound 
system is revisited in various regimes. First, weakly gravitating Newtonian systems are considered, 
followed by various exact solutions describing a relativistic central object embedded in a Friedmann 
universe. It is shown that the "all or nothing" behaviour recently discovered (i.e., weakly coupled 
systems are comoving while strongly coupled ones resist the cosmic expansion) is limited to the de 
Sitter background. New exact solutions are presented which describe black holes perfectly comoving 
with a generic Friedmann universe. The possibility of violating cosmic censorship for a black hole 
approaching the Big Rip is also discussed. 



INTRODUCTION 

The issue of whether a planet, a star, or a galaxy ex- 
pands following the rest of the universe is a problem of 
principle in general relativity that still awaits a defini- 
tive answer. The effect of the cosmological expansion on 
local systems such as the Solar System has a long his- 
tory dating back to the 1933 paper by McVittie [l] intro- 
ducing a spacetime metric that represents a point mass 
embedded in a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe. Later work by Einstein and Straus 
[21 introduced the Swiss-cheese model which is, however, 
unable to describe the Solar System 0, Q] and is unlikely 
to be extended to non-spherical systems [a, lal. Many 

1, El, m 



21|, 122, 123|, 1241 presented 



papers in the followin g ye ars 

m, Q E [S [13, [3, 111, H 

contradictory results, casting a shadow of ambiguity on 
the problem (see [3,[2^ for brief reviews). The most pop- 
ular model consists of a test particle in a quasi-circular 
orbit around a Newtonian central object. It appears that 
many of the contradictory results are simply due to the 
use of different or unphysical coordinates [l3|. Moreover, 
a quantitative answer to the problem of how much the 
cosmic expansion affects local dynamics differs accord- 
ing to the type of local system considered. If the FLRW 
metric is an adequate model of spacetime down to small 
scales, then weakly gravitating systems of size small in 
comparison to the Hubble radius Hq^ do participate in 
the expansion, but the effect is so small to be completely 
negligible for practical purposes. When the size of the 
weakly gravitating system becomes a larger fraction of 
the Hubble radius, the cosmic expansion plays a signif- 
icant role in the dynamics; this is the situation of large 
scale structures [1, [13, |23 ■ 

A recent paper by Price [22| studies a classical atom 
in a de Sitter background, with arbitrary strength of the 
coupling between an electron and the central charge. A 
new result is the "all or nothing" behaviour: if the cou- 
pling is weak the "atom" is comoving with the rest of the 
universe, while if the coupling is very strong the "atom" 



is only slightly perturbed by a transient and does not 
expand [22]. This work breaks free of the standard as- 
sumption of previous literature that the coupling (of a 
gravitationally, instead of electrically, bound system) is 
weak. However, it has two fundamental limitations: first, 
the cosmological background is restricted to be de Sitter 
space, which is very special: in fact, the de Sitter met- 
ric can be put in static form, which may explain why 
strongly bound systems in this background do not ex- 
pand. Second, the classical atom is not an adequate 
model when the local energy density and stresses of the 
central charge grow and induce local deviations from the 
cosmological metric. In these situations exact solutions 
of the Einstein equations are needed to describe both the 
relativistic central object with a strong local field and the 
surrounding universe. Solutions of this kind could also be 
useful in studying the evolution of primordial black holes 
[i^t regarded as probes of the early universe 27, |2^. 
New and existing exact solutions of this kind are studied 
in Sees. 3 and 4. 

There are also more modern and perhaps more com- 
pelling motivations to study the effect of the cosmological 
expansion on local physics. It is now well known from the 
observations of supernovae of type la that the expansion 
of the universe is accelerated [29| . Marginal evidence for 
an equation of state parameter w = ^ < —1 (where p 
and P are the energy density and pressure of the cos- 
mic fluid, respectively) has led theorists to take seriously 
into account the possibility of a Big Rip singularity at 
a finite time in the future (30| . Various authors have 
studied how local systems (clusters, galaxies, stars, etc.) 
are teared apart as the Big Rip is approached 21 , [s^ . 



In this situation the catastrophic cosmological expansion 
is not merely a perturbation of the local dynamics, but 
dominates it. 

The current inability to explain away the Pioneer 
anomaly has led some authors to attribute it to the ef- 
fect of the cosmological expansion, although this possi- 
bility seems to be ruled out [20]. Finally, independent 
motivation for studying the interplay between local and 
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cosmological dynamics comes from another problem of 
principle in general relativity. If a universe dominated 
by phantom dark energy, which violates all the energy 
conditions and causes w = P/p to he less than —1, is 
heading toward the Bi g R ip and all bound systems are 
gradually ripped apart [30| , it is legitimate to ask what is 
the fate of the most strongly bound local object, namely 
a black hole. Does the strong local field resist the expan- 
sion (as suggested by extrapolating Price's work [ISDi or 
does the horizon expand and disappear, exposing the cen- 
tral singularity before the Big Rip is reached? This would 
entail violation of cosmic censorship in its cosmological 
formulation 'sT] and would constitute a (further) argu- 
ment against phantom energy. An answer comes from 
Ref. [32| , in which accretion of a phantom test fluid onto 
a Schwarzschild black hole is studied. The results are 
extrapolated to a gravitating fluid and the conclusion is 
reached that the black hole decreases its mass due to the 
fact that the gravitating energy density accreted P + p 
is negative. Accretion proceeds until the horizon disap- 
pears together with the central singularity before the Big 
Rip is reached (s^ . Although the extrapolation from a 
test to a gravitating fluid is quite plausible, it would be 
preferable to base the conclusion on an exact solution 
of the Einstein equations displaying accretion of cosmic 
fluid. A step in this direction is taken with new solutions 
presented here. 

In this paper we examine various exact solutions repre- 
senting strong field objects in a cosmological background: 
they include the McVittie metric the Schwarzschild- 
de Sitter black hole, the Nolan interior solution [s^l, a so- 
lution found recently by Sultana and Dyer [l9l |. and new 
exact solutions that are perfectly comoving. It turns out 
that the "all or nothing" behaviour discovered by Price 
[2^ persists in the Schwarzschild-de Sitter black hole but 
it is a peculiarity of the de Sitter background adopted 
and more general FLRW backgrounds do not allow for 
it. Participation of a local object, even strongly bound, 
to the cosmological expansion seems to be the general 
rule, a conclusion supported by various exact solutions. 

The plan of this paper is as follows: Sec. 2 studies 
a Newtonian quasi-circular orbit and the effect of the 
cosmic expansion upon it, making clear the peculiarity 
of de Sitter space even for this kind of systems. Sec. 3 
examines known and new exact solutions, while Sec. 4 
contains a discussion and the conclusions. We adopt the 
notations of Ref. [3J| . 



A NEWTONIAN OBJECT EMBEDDED IN A 
FLRW UNIVERSE: QUASI-CIRCULAR ORBITS 

In the literature, the most common line of approach 
to the problem of the effect of the cosmic expansion on 
local systems is to consider a spherical Newtonian ob- 
ject of mass M and a test particle in a circular orbit of 



radius r around it, and then "switch on" the cosmolog- 
ical dynamics as a small perturbation of this two-body 
problem (elliptical orbits were considered in Ref. 
Motivated by the current model of our universe (and by 
simplicity), and following most authors, we consider as 
the background a spatially flat FLRW universe described 
by the line element 



(1) 



in comoving coordinates. The evolution equation for the 
physical radial coordinate of the otherwise circular orbit 
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r = — r 

a r 



GM 



(2) 



where L is the (constant) angular momentum per unit 
mass of the test particle and an overdot denotes differen- 
tiation with respect to comoving time. This equation is 
derived in several ways in the many papers on this sub- 
ject: they range from heuristic derivations (e.g., [tL \2^) 
to calculations using the geodesic deviation equation in 
a locall y in ertial frame of the cosmological metric ^ 
Til . 24, [3^. This derivation uses Fermi normal coordi- 
nates, regarded as the physical coordinates connected to 
a freely falling observer in the cosmological gravitational 
field. Other derivations of eq. ^ [sj, l2l| use various ap- 
proximations to equations for timelike geodesies in the 
McVittie metric in the limit in which a central ob- 
ject produces only a small deviation from the cosmolog- 
ical background. If the current era in the history of the 
universe is considered, the term dr/a on the right hand 
side of eq. (I2|) is a small perturbation when r << Hq^, 
where Hq is the present value of the Hubble parameter 
H = d/a. This correction becomes increasingly impor- 
tant as the size of the orbit increases in comparison with 
the Hubble radius H^^ , for example going to galaxy clus- 
ters and superclusters, for which it is certainly not neg- 
ligible BSEIii. 

Note that, for a decelerating universe, the term - r is 
negative and is considered as such in pre-1998 literature, 
while it gives a positive contribution to f in an acceler- 
ated universe, which is the case of the present epoch, as 
is now well known from the study of supernovae of type 
la at high redshift [2^. 

Whether dr/a can be treated as a perturbation or not, 
both the central object and the FLRW background are 
assumed to be spherically symmetric, which yields con- 
servation of the angular momentum per unit mass L of 
the test particle, 



(3) 



in spherical coordinates. This equation holds true in later 
sections in which we consider exact solutions of the Ein- 
stein equations describing a strongly gravitating, spher- 
ically symmetric, central object embedded in a FLRW 
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universe. It follows that, if r{t) increases with time, the 
angular velocity (p{t) decreases. The test particle has to 
cover a larger linear distance to attempt to close its orbit, 
which would be circular in the absence of the cosmologi- 
cal perturbation. 

A different point of view is the one, found in recent 
literature, in which it is assumed that the universe is 
dominated by phantom energy with equation of state pa- 
rameter w = P/ p < —1 [3^, and is heading toward a Big 
Rip in which the scale factor a(t) diverges at a finite time 

A long time before the Big Rip, the term ^ r can be 
treated as a small perturbation. In flat space, at the 
Newtonian level, this term is absent, Kepler's third law 

^2 ™3 — 



Then <p{t) (po a.st ■ 



rip ' 



as the Big Rip is approached 



and r grows without bound (but comoving) , the angular 
motion slows down and freezes. Strictly speaking, the or- 
bit is never "disrupted" before the spacelike singularity 
is reached; it just participates in the cosmic expansion 
that is accelerating catastrophically. In this sense, it is 
not true that bound systems become unbound: this is a 
rather misleading sentence often echoed by the media. If 
one wants to insist on the use of this terminology, the 
meaning of "bound" and "unbound" system should be 
clearly defined. For example, one may think of (arbi- 
trarily) setting the threshold between "bound" and "un- 
bound" when, in eq. the cosmological term dr/a be- 
comes of the order of the other two terms —GM/r^ and 



yields ip'^r'^ — GM, and the two terms —GM/r^ and I? jr". The condition drja GM/r^ can be expressed 



I? jr^ in eq. ^ cancel each other, leading to circular 
orbits r ^constant. When the perturbation f t" is in- 
troduced, this is no longer true and this perturbation 
changes the (otherwise circular) orbit. When the Big 
Rip is approached as t ^ i~p and the central object is 
weakly gravitating (Newtonian) , the term - r dominates 
over the terms in r^^ and r~'^, leaving 



(4) 



as the asymptotic evolution equation for the physical ra- 
dius of the orbit. The solution of eq. is r oc a(i) or, the 
orbit becomes comoving with the cosmic substratum. For 
simplicity, assume that the phantom energy dominating 
the cosmic dynamics and causing the Big Rip has con- 
stant equation of state P = wp^ where w < — 1. Then, 
the scale factor is 



by saying that the time scale t ~ is of the order of 
the free fall time scale of a fictitious region of radius r 



containing the mass M, i.e. 



where f £ is the 
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escape velocity. Or, in other words, the energy density 
~ M/r^ of this fictitious region equals the cosmologi- 
cal density ^ d/a which grows in a phantom-dominated 
universe. A more precise characterization of when a 
"bound" system is "disrupted" is given by Nesseris and 
Perivolaropoulos 121]. These authors consider the effec- 
tive potential V{t, r) of the one-dimensional equation of 
motion of the test particle and determine when its min- 
imum disappears by solving numerically the equation 
dV/dr = (the location of this minimum depends on 
time). This procedure corrects our order of magnitude 
estimate by a factor ~ 3. The effective potential used is 
the subject of the next subsection. 



a{t) = ao (trip - t)" 



3(u; + l) 



< 



and eq. (|4]) reduces to 



a (a — 1) 



(5) 



(6) 



(7) 



{irip ^) 

which has the general solution 

r{t) = A {t„p -tf+B {t„p - t)'"" , 

where A and B are integration constants. Since 1 — a — 
^1^^}^ > 0, the second term on the right hand side of 
eq. ([7]) becomes negligible with respect to the first one 
as the Big Rip is approached. Therefore, the solution 
d?]) reduces to r{t) cx a{t) as i — > i"^, i.e., the putative 
circular orbit becomes comoving. 

The angular motion is obtained by integrating eq. ([31), 
which yields 



^^L^3L(«. + 1) 



1 - 3w 
1) L 



3(w . - /, 

I — 6W A"^ 



A^ (i„p-t)3(»+i) +AB 



(8) 



Effective potential, Lagrangian, and Hamiltonian 



By rewriting eq. ^ as 



GM 



^3 



dV 
dr 



(9) 



integrating with respect to r, and setting to zero an ar- 
bitrary integration function of time, one obtains the ef- 
fective potential 



a 2 GM 



2r2 



7^ ■ (10) 



In a general FLRW spacetime this effective potential for 
the bound system test particle-central object depends 
on time because of the cosmological term —^r^ and 
one can say that this system exchanges energy with the 
cosmological background. However, in the special case 
of a de Sitter background described by the scale factor 
a(t) — aoexp(_ffot) with constant Hq, this term is time- 
independent and reduces to — — no energy is 
then exchanged between the two-body system and the 
cosmological background, and the energy of the former 
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is constant. This simplification was noted by Price [22l |. 
who restricted his study of the effect of cosmological dy- 
namics on local systems to a de Sitter background. 

Given the effective potential PU)) . it is straightforward 
to deduce the effective Lagrangian 



GM 



2r2 



2a 



(11) 



which reproduces eq. ^ through the Euler-Lagrange 
equation 



d_ 

di 



dC 

dr 



dC 

dr 



and the effective Hamiltonian 

GM 







2r2 



2a 



(12) 



(13) 



where Pr = dC/dr = r and 

. dn 

OPr 



dn 

dr 



(14) 



This Hamiltonian is, of course, time-dependent for any 
FLRW background that is not a de Sitter space: 



dn _ _dC _ dV 
'dt ~ ^~dt ~ 'dt 



(15) 



For a de Sitter background a = aoe ° with Hq — v A/3 
(where A > is the cosmological constant), the effective 
potential (fTO| is consistent with the weak-field limit of 
the Schwarzschild-de Sitter (or Kottler) metric, given by 
the line element 



1 - 

r 



2GM Ar 



dr^ 



(16) 



in static coordinates, where dQ^ = d0^ + sin^ 9 dip^ is respect to the Hubble radius H^^ ^. When radii r ^ Hq 



the line element on the unit 2-sphere. The analogue of 
the Newtonian potential $jv for a test particle can be 
read off the (0, 0) component of the Kottler metric goo — 
- [1 + 2<i>Ar(r)], which yields 



$Ar(r) 



GM 



Ar2 
~6~ 



(17) 



On the other hand, the cosmological part of the potential 
PH)) in the equation of motion ^ of the test particle in 
the de Sitter background with Hq = -\/A/3 is 



Hir^ 



2a 



(18) 



_Ar2 

This is not a coinci- 



which is consistent with eq. ([T7| . 
dence because eq. ([2|) can be derived as a special approx- 
imation of the timelike geodesic equation of the McVittie 
metric Jj], which reduces to the Kottler metric for a de 
Sitter background (see Sec. 3.3). 

By adding the centrifugal potential term ^ one re- 



covers the effective potential ([TOl) for the equivalent one- 
dimensional problem and V{r) — <i>jv(r). 

Orbits of constant radial coordinate 

It is time to comment on the physical meaning of the 
radial coordinate r. As shown in Ref. [l^l, this coincides 
with the proper radius when rHo << 1 and receives cor- 
rections of higher order when r is larger and larger with 



are considered, r assumes the meaning of comoving radial 
coordinate in the FLRW metric (HI). However, when the 
field of the central object is strong, neither of the above 
describes precisely the meaning of r for r << Hq^ . 

The question of whether orbits of constant radial co- 
ordinate r exist is of some interest. Orbits of constant 
radius are found by Bonnor [3] who, similarly to Price 
j22,] , considers a classical atom embedded in a cosmolog- 
ical background. Two situations can be distinguished, 
i) The cosmological background is not a de Sitter 
space. Then, d/a depends on time and, imposing r = 
ro =constant, one obtains 



GM L2 



(19) 



The first term on the left hand side depends explicitly on 
time, while the remaining terms are time-independent, 
hence this equation can not be satisfied and orbits of 
constant r do not exist in this case, 
ii) The cosmological background is de Sitter space. Then 
imposing r — ro =constant yields the quartic equation 
for ro 



HM - GMro + 







(20) 



This equation can, in principle, have real solutions un- 
der conditions which are discussed in the next subsection. 
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Test particle around a Newtonian central object in a 
de Sitter background: phase plane analysis 

The equation of motion of an electron in a classical 
atom embedded in a de Sitter background is analogous to 
our eq. ^ and is solved numerically by Price [1^. Con- 
trary to previous authors, Price does not restrict himself 
to considering weak couplings of the electron, which is 
the equivalent of the situation considered so far in our pa- 
per, of a test particle in the field of a Newtonian, weakly 
gravitating central object embedded in a cosmological 
background. Price considers instead arbitrarily strong 
coupling of the electron to the central charge and dis- 
covers an "all or nothing" behaviour: if the coupling is 
weak the electron trajectory becomes comoving with the 
de Sitter substratum while, if the coupling is strong, the 
evolution of the orbit exhibits a transient after which it 
is essentially unperturbed. A critical value of the an- 
gular momentum separates these two behaviours. This 
"all or nothing" feature went undetected in the (abun- 
dant) previous literature which did not break free of the 
weak coupling assumption. In our formalism, allowing 
for an arbitrarily strong gravitational field due to the 
central object means making this object relativistic and 
one must leave the regime in which the latter is merely 
a perturbation of a cosmic substratum, and move to a 
fully relativistic regime by studying exact solutions de- 
scribing a strongly gravitating central object embedded 
in a FLRW background. This will be done in the next 
section. Here we limit ourselves to weak coupling and we 
provide a phase space analysis of the equation of motion 
analyzed numerically by Price 22], in the analogous sit- 
uation of a two-body system completely ruled by gravity. 
This phase space analysis is still missing in the literature. 

By introducing variables x = r and y = r (not to 
be confused with Cartesian coordinates in the de Sitter 
background ([1])), eq. ([2]) is written as the autonomous 
dynamical system 



X = y 



y = H^x 



GM 



x^ 



(21) 



(22) 



for X > and any real value of y. We assume L 7^ 
(the case L = will be discussed later). Equilibrium 
points, if they exist, correspond to orbits of constant ra- 
dius (a;, y) = (ro, 0) of the kind discovered by Bonnor 
The search for these fixed points is equivalent to solving 
eq. (I2ni), or 



i!{x) = Hlx^ - GMx + = 



(23) 



The function ip{x) is represented by a quartic parabola 
with its concavity facing upward, therefore there can be 
solutions of eq. (I23p only if the minimum ^min of ipi^) is 
non-positive. If ipmin = there are two coincident roots. 



while if 4'min < there are two real distinct roots, and 
no real roots exist if ■(/'mm > 0. The study of the first 
derivative dtp/dx establishes that ip{x) has the absolute 
minimum 



-0 {Xr, 



2 3z 
-"0 



at 



where 



/GMHo\''^^ 



(24) 



(25) 



(26) 



is a dimensionless variable. It follows that V'mm > 0, 
hence there are no orbits of constant r if 



L>V3 



GM 



2/3 



(27) 



There is a single orbit of constant r = Xmin if ^ = L^, 
and there are two orbits of constant radii ri 2 with ri < 
^min < ^2 whcn L < Lc- Thcsc are all the fixed points 
(xqjO) of the dynamical system (pT|) and ((22)) . In order 
to assess the stability of these fixed points let us consider 
perturbations described by x{t) = xq + dx{t), y{t) = 
5y{t). Eqs. ([2T|) and ([22]) yield the evolution equation for 
the orbital radius 



5f 



0. 



= -r - 3i/, 







(28) 



Linear stability corresponds to up' > 0. Orbits of con- 
stant r exist when L < Lc, with the equality correspond- 
ing to a single orbit of radius tq and the strict inequality 
to two orbits of radii ri^2 with ri < Xmin < 7*2 • The 
stability condition of an orbit of constant radius is 
equivalent tor* < . / J" . In the case of a single orbit 

Y v3 Ho 

it is L = Lc and = rg = Xmin < Xmin and the stability 
condition is satisfied. In the case of two distinct orbits 
it is ri < Xmin < ^2 and therefore the outer orbit is un- 
stable. The inner orbit is stable only if ri < ^/ J"^^ . 

V V 3 -"0 

Since (. = = If - GAf < when 

L > Lc and ipi^) is a decreasing function between x — 
and Xmin, while i^i^i) — 0; it must be ri < ^ and 
therefore the inner orbit is stable. 

At large values of x, eq. ([22)1 reduces to the asymptotic 
equation y ~ H^x, which has the solution 



{x{t),y{t)) = {x,e"°\Hox,e"°') 



(29) 



represented by the line y = H^x in the {x, y) plane. This 
solution is an attractor for large values of x. In fact, by 
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perturbing the solution as described by x{t) = x^e^"* + 
Sx{t) and y{t) = iJox^e^"' + Sy{t) one obtains, to first 
order in the perturbations, 6x = Hq5x+ ... with solution 
6x{t) = SqC^"* (with 60 a constant). Therefore, the ratios 



dx{t) 
a;*e-f^"* 



xo 



JT^ - - (30) 



«1) 



stay small if they start small ( 

The qualitative picture of the phase space can be com- 
pleted as follows. The system ((2T|) and (|22l) has the first 
integral 



E : 



2 



H^x^ 



GM 



(31) 



corresponding to the energy of the test particle, which is 
conserved in a de Sitter background. 

• li L > Lc there are no fixed points and y = Hqx is 
an attractor for large values of x. Then y — ip{x)/x^ > 
(the equation -0 = has no real roots for L > Lc). 
Hence y = x always increases and either y +00 or y{t) 
has a horizontal asymptote with y{t) —> yo sls t ^ +00. 
If y = i — > +00, then either x(t) +00 or x{t) has a 
vertical asymptote with x xq and x — s- +00. There 
are, in principle, four possibilities. 

1) x{t) — !■ +00 and y{t) +00; then the orbit of the 
solution is necessarily captured in the attraction basin 
of the attractor y — Hqx. 

2) x{t) — > +00 and y{t) yo; this is not possible 
because as a; +00 it is y = Hqx — )■ +00. 

3) x(t) — > a;o and y{t) —>■ +00 {x(t) has a vertical 
asymptote x +00); then the energy E given by 
eq. ([5T|) diverges at late times, which is not possible 
because E is constant. 

4) x(t) — > xq and y(t) — )■ yo; this means that there is an 
attractor point, then y — H^x — GMx^^ + L?x^^ 
as t ^ +00 then must be a root of the equation 
y — Q but this is not possible because this equation has 
no real roots when L > Lc. This case is not possible. In 
summary, for L — Lc all the trajectories of the solutions 
in phase space go to the attractor y = Hqx as t —^ +00, 
i.e., all physical orbits become comoving. 



• If L = Lc there are the attractor y 

and the fixed point xq = Xmin = (^^^ 
is the only real root of the equation y 



Hqx as a: — > +00 

1/3 

,2/0 = which 
: 0. Since there is 



only one such root, a priori either y > or y < along 
any orbit that does not coincide with the fixed point, 
without the possibility of changing sign during the evolu- 
tion. However, we know that y ~ ip(x)/x^ > \/x ^ xq 
then y = X is always increasing for x ^ xq. By repeating 
the reasoning of the case L > Lc, we have now two 
possibilities: 

1) x(i) +00 and y{t) +00, in which case the orbit 



of the solution in phase space gets captured by the 
attractor y = Hqx at infinity. 

2) ix{t),y{t)) ^ (xo,0) with y{t) Q- so y = x < 
and x{t) is always decreasing to Xq. Then, a; > 
with i < 0; the point (a;o,0) is an attractor point. 
Since y > 0\f{x,y) ^ (xo,yo) there are no periodic 
orbits. The cases in which {x{t),y{t)) {+00, yo) 
or {x{t),y{t)) (xo,-l-oo) are excluded as previously 
discussed for L > Lc. In summary, for L = Lc the orbits 
of the solutions either go to the attractor point or to the 
attractor at infinity y = Hqx. 

• If < L < Lc there are the equilibrium points 
(^1.2:0) and the attractor at infinity y = Hqx with 
ri < Xmin < 7*2 • There are two real distinct roots of the 
equation y = ip(x)/x^ = and it is y > for x < ri and 
for x > r2, while it is y < for ri < x <r2. 

For X < ri the function yit) = i is increasing and 
either y ^ +00 or y(t) — > y^ (horizontal asymptote of 
y(t)). If y ^ +00 it must be because a; ^ 0, the familiar 
situation in which the Newtonian centrifugal potential 
L^x~^ dominates and repels the particle to larger x. If 
y(t) has a horizontal asymptote, y — > yo, then y ^ and 
the attractor point (ri , 0) is approached. 
For ri < a; < r2 it is y = tp{x)/x^ < so y{t) is always 
decreasing and either y(t) —^ —00 or y(t) has a horizon- 
tal asymptote y —t y^ . The first situation is not possible 
because it would imply that the energy E diverges, while 
it is instead forced to be constant. Therefore, y{t) must 
have a horizontal asymptote and y ^ 0. Then the at- 
tractor point (zero of the function ili{x)) is approached. 
For a; > r2 it is y > and either y{t) -\-oo or y{t) — > y^ 
(horizontal asymptote). If y ^ -l-oo there are, in princi- 
ple, two possibilities: 

1) x{t) +00 and y{t) +00; then the orbits of the so- 
lutions are captured by the attraction basin of y = Hqx. 

2) x{t) stays finite and y{t) +00. Again, this would 
give E +00, which is excluded. 



Zero angular momentum 

Finally, let us consider the special case L = in which 
the radial coordinate of the test particle satisfies 



r = Hir - 



GM 



(32) 



Only one "orbit" of constant physical radius ro = 

1 /3 

{GM/H^) ' exists in this case. This situation is not 
possible in flat space and it corresponds to the cosmo- 
logical "force" H^r exactly compensating the attractive 
force of the central object —GM/r^. As is clear from 
the physical point of view, this position of equilibrium is 
unstable because an arbitrarily small radial displacement 
will move the test particle to a region where one of the 
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two forces is dominant. The value tq — {GM/Hq 
corresponds to a maximum of the effective potential 



1/3 



V{r) = - 



GM 



(33) 



In fact, setting dV/dr 
and d^V/dr' 



= -3H^ < 



~ reproduces this value of r 
This simple solution is 



apparently missed in previous literature. 



Purely radial motion in a general FLRW background 



simplified and perhaps is not the best physical system in 
which to study the competition between local attraction 
and cosmological expansion. The next step would be to 
consider the expansion of a Newtonian star embedded in 
a FLRW universe. We skip this step and, in the next 
section, we consider instead a relativistic star embedded 
in a FLRW universe, and its Newtonian limit. 



A STRONGLY GRAVITATING OBJECT IN A 
FLRW BACKGROUND 



One can generalize to any FLRW background the 
search for solutions with purely radial motion satisfying 



a GM 
— r TT- 



(34) 



Orbits of constant r are impossible if the background is 
not de Sitter space. We consider universes with scale fac- 
tor given by a power-law a(t) = oq (with p a constant). 
Then, it is easy to find the power-law solution 



r*it) 



GM 



-,1/3 



p{p-l) + 



.2/3 



2/3 



(35) 



This solution exists only for values of the constant p satis- 
fying p < 1/3 or p > 2/3, which follows from the require- 
ment that the denominator of the constant Vc in psp be 
positive. The solution ([55)1 is unstable with respect to 
small radial perturbations, In fact, if r{t) = r^,{t) + Sr{t), 
using eq. (j35p one finds the evolution equation for the 
linear perturbations Sr 



Sr 



3p(p-l) 



(36) 



which has solutions Sr{t) = (Jq with (Jg a constant and 



1± 




(37) 



In the allowed range of values of p, the discriminant A = 
12p(|5- 1) + ^ satisfies VA > 1/3, which yields /3+ = 

^1 + \/A^ /2 > 2/3 and the corresponding mode 5r+(t) 
satisfies 



fr+(0 
r*(t) 



(38) 



Because of the positive exponent /3+ — 2/3, this ratio 
grows when t +oo and therefore the solution ([35)1 is 
unstable. 

To conclude, one recognizes that the de Sitter universe, 
although simpler to study than a general FLRW back- 
ground, is a very special case. Furthermore, the simple 
problem of a test particle in a circular orbit is extremely 



The study of a classical atom in a de Sitter background 
by Price [22| shows that systems that are strongly bound 
"resist" the cosmological expansion and are only per- 
turbed a little by it. This leads one to believe that this 
situation would carry over to the case of a gravitation- 
ally bound system and that a strongly gravitating central 
object, such as a black hole, would not be perturbed at 
all by the cosmic expansion. Whether this is true or not 
is the first question addressed in this section. Second, 
one would like to know whether this time-independent 
behaviour of a strongly gravitating central object (if it 
really carries over from Price's classical atom to the grav- 
itating system) is peculiar of de Sitter space, or if it is 
valid in general FLRW backgrounds. After all, de Sitter 
space (other than the trivial Minkowski space with a = 1) 
is the closest to a static space, and (a portion of) it can be 
expressed in static coordinates; therefore, it could be too 
special to derive general conclusions. Third, it is claimed 
in the literature that in a universe approaching the Big 
Rip all bound objects (galaxies, stars, atoms, etc?) are 
ripped apart before the singularity is reached [30[ . What 
about a black hole horizon? If the expansion to a Big Rip 
tears apart a black hole horizon and a naked singularity 
appears, the cosmological version of cosmic censorship 
|3ll | is violated, the implications for black hole thermo- 
dynamics are non-trivial, and the phantom energy caus- 
ing the Big Rip would be questioned further from the 
point of view of fundamental physics. Regarding this 
problem, it is claimed that a black hole accreting phan- 
tom energy with P < —p in a universe approaching the 
Big Rip decreases its mass and the horizon disappears 
together with the central singularity before the Big Rip 
is reached [s^] ■ This phenomenon avoids the violation of 
cosmic censorship by eliminating the central singularity 
altogether. New exact solutions can help clarifying this 
issue. 

In all the situations mentioned above one can heuristi- 
cally see the "all or nothing behaviour" as resulting from 
the competition between two strong fields, the cosmo- 
logical one and the local one due to the central object. 
It is conceivable that one of the two is locally stronger 
than the other and "wins" but, on the other hand, many 
examples of black holes are known which, placed in a 
strong external gravitational, electric, or magnetic field, 



8 



have their horizons stretched 11 



H, H S, EE 53, S, H 0, E 



m, [si, EE El, Si 



52| 



In this section we want to give a more precise mean- 
ing to these naive considerations, and a general picture 
of these phenomena that is not limited to a de Sitter 
background space. Exact solutions of the Einstein equa- 
tions representing strongly gravitating objects embedded 
in a FLRW universe are needed. We neglect semiclassical 
Hawking radiation in the following. 



The McVittie solution 

An obvious starting point is the McVittie metric intro- 
duced in 1933 with the explicit purpose of investigatijig 
the effects of the cosmic expansion on local systems l] . 
The line element is 



1 - 



m{t) 
2f 



m(t) 
2r 



■dt^+a\t) 1 + 



m{t) 



(39) 



in isotropic coordinates, where the function m,[t) of the 
comoving time satisfies the equation [l[ 



m 
m 



a 
a 



which yields 



m{t) 



a{t) 



(40) 



(41) 



where niH is a constant representing the Hawking- 
Hayward quasi-local mass [HE HJ] (see Refs. [1, [THj for a 
discussion) . Eq. ([40)l was imposed by McVittie to explic- 
itly forbid the accretion of cosmic fluid onto the central 
object (assumption e) of (H). It corresponds to Gj = 0, 
which in turn implies that the stress-energy tensor com- 
ponent Tq — and there is no radial flow. In modern 
language, eq. (|40)) corresponds to the constancy of the 
Hawking-Hayward mass, ttih = 0. It is important to rec- 
ognize TTiH as the physically relevant mass (eventually re- 
lated to the physical size of the horizon or of the central 
object) in order to avoid making coordinate-dependent 
statements on the mass and size (cf., e.g., Refs. [l8l[55|). 
or temperature [1^ of the central object. m{t) is just a 
metric coefficient in a particular coordinate system. 

While there is little doubt that the McVittie metric 
represents some kind of strongly gravitating central ob- 
ject, its physical interpretation is not completely clear 



and is still debated today IE, 12, 15, [lE|. This metric 
reduces to the Schwarzschild solution in isotropic coor- 
dinates when a = 1 and to the FLRW metric if m = 0. 
However, in general, the metric (|39p can not be inter- 
preted as describing a black hole embedded in a FLRW 
universe because it is singular on the 2-sphere f — m/2 
(which reduces to the Schwarzschild horizon if a = 1) 
lE IE iE| and this singularity is spacehke 15, lE|. It is 
claimed that the McVittie metric describes a point mass 
located at f = (which is another, expected, singularity 
of the metric) embedded in a FLRW universe. However, 
this point mass is, in general, surrounded by the singu- 
larity at f = m/2, which is difficult to interpret. This 
singularity was studied in Refs. [3, [13, [H]. Nolan i 
showed that this is a weak singularity in the sense that 
the volume of an object falling onto the f = m/2 surface 
is not crushed to zero, and therefore the energy density 
of the surrounding fluid is finite. However, the pressure 



P = - 



1 



2H{l 



1 m. 

^ 2r 



(42) 



diverges at r = m/2 together with the Ricci scalar 
R = 87rG(3P-p) [IE [3 [IE [11. Notwithstanding 
this, there is a situation in which the singularity dis- 
appears and the surface f = m/2 describes a true black 
hole horizon: this happens when the background FLRW 
is de Sitter space. Then H = and the second term on 
the right hand side of the expression (j42|l — the only one 
causing P to diverge — is absent (this point was also noted 
in Ref. [lEl)- A possible reason for the disappearance of 
the singular surface in the de Sitter case is discussed be- 
low. Similar problems affect the charged McVittie met- 
ric 18] and the solutions of Thakurta [13], Vaidya [sst - 
Patel and Trivedi [3] representing rotating black hole- 
like objects in a cosmological background ([551 — see also 
Refs. [4, 56]). 



The Nolan interior solution: a relativistic star in a 
FLRW universe and its Newtonian limit 

It is of interest to study the behaviour of a relativis- 
tic star embedded in a FLRW background with respect 
to the problem of local physics versus cosmological ex- 
pansion. The Nolan interior solution [s^] describes a rel- 
ativistic star of uniform density in such a background. 
The metric is 



9 



^ rn 

^0 



1 7n 

4fo 



(43) 



in isotropic coordinates, where rg is the star radius, 
— = — - (the condition forbidding accretion onto the star 
surface), and < f < tq. The interior metric is regular at 
the centre and is matched to the exterior McVittie metric 
at f = fg by imposing the Darmois-Israel junction condi- 
tions. The energy density is uniform and discontinuous 
at the surface f = tq, while the pressure is continuous. 
These quantities are given by [s^l 



Pit) 



P{t,f) 



8ttG 



8ttG 



3H^ 



6m 



-3H^ - 2H 



(44) 



4ro I r, 



3m" 



m \ mr^ 

4fn 



(45) 



The Nolan interior solution is a special member of Kus- 
taanheimo's family of shear-free solutions 59] that gener- 
alizes the Schwarzschild interior solution of uniform con- 
stant density to the case of a time-dependent cosmolog- 
ical background. By setting a = 1 one recovers familiar 
expressions for the Schwarzschild interior solution 34j . 



The energy density is always positive and the condition 
P > imposed by Nolan coincides with d + 3a^/2 < 

Let So(i) = {{t,r,6,ip) : r = r^} be the surface of 
the star at time t; by construction, the metric on this 
2-sphere coincides with the McVittie metric 



1 - 



m(f) 
2fo 



m{i) 
2fo 



■dt^ + a\t) 1 + 



m{t) 
'2f7 



(46) 



The proper area of Eq is 



where gab 



(So) 



m(t) 

; ^ ^ (47) 

is the metric on Sg at a fixed time t and 



(/So is its determinant. By using the Schwarzschild cur- 



vature coordinate r = r (l 



^) , one has 



^So(0 



47ra2(t)r^. (48) 
The star surface is comoving with the cosmic substratum 
and the proper curvature radius of the star is rphys (t) = 

a{t) fg ^1 -I- . Therefore, we have a local relativistic 
object with strong field which is perfectly comoving at 
all times: in this case the cosmic expansion "wins" over 
the local dynamics. 

It is interesting to compute the generalized Tolman- 
Oppenheimer-Volkoff equation fs^l valid for this crude 
star model and to derive the first order correction to the 
Newtonian equation of hydrostatic equilibrium. For a 
perfect fluid described by the stress-energy tensor Tab — 
[P + p) UaUb + Pgabi the covariant conservation equation 
V^Tab ~ splits into the two equations 34 1 



u^V^p^ {P + p) = 



h\d,P+{P + p) u'^W.Ua^Q 



(49) 



(50) 



where u° is the fluid four-velocity and hab = Qab + UaUb 
defines the projector ha'^ onto the 3-space orthogonal to 
In comoving coordinates it is u'^ oc 5'^'^ and the nor- 
malization u'^Uc = — 1 yields 



or u = |goo 



1 



2fo 



)( 



1 



fa 



1 



(1,0,0,0) , (51) 



4ro J rg 

Eqs. (gSl) and dMl) then yield 



dp 
dt 



3H{P + p){l 



ro 



2 1 



2fo J 



2fl 



2rl 



0, 



(52) 



which generalizes the well-known conservation equation 
p + 3H {P + p) =0 valid in a FLRW universe, to which 



it reduces in the limit m — > 0. For the Schwarzschild 
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interior solution there is no equation analogous to (1521) 
because = for the static Minkowski background and 
p is static. 

Eq. ([50)1 can be rewritten as 



(53) 



dP , „ X mr 
7P+(P + P) — 
or r}] 



1 



m \ 
2ro J 



1 



Iff 



1 



1 



4fo i 



(54) 

In the Newtonian limit m/r, m/ri) « 1, P << p, r r, 
this equation reduces to 



dP 

dr 



d<i> 



N 



dr 



where p — m rg) 



and — 







2fA 



dP d^N 3 . , , ^ , 







(56) 



A= dOdip^ , 



(57) 



where gab Is is the restriction of the metric tensor to this 
2-sphere and g^. is its determinant. Eq. p9|) yields 



or, upon use of eq. (j41|) . 

^ = 16n rnjj 



(58) 



(59) 



The Schwarzschild radial coordinate corresponding to the 
horizon is the curvature coordinate 



By setting the index c — 1 and computing the covariant 
derivative one obtains 



= 



(55) 



is the Newtonian 



potential. This expression for the density can also be ob- 
tained from eq. (|44ll by setting a = 1 and using the cur- 
vature radius. The first order correction to the equation 
of hydrostatic equilibrium for a spherically symmetric, 
uniform density star, is given by 



The Schwarzschild-de Sitter black hole 

When the de Sitter space is chosen as the background, 
it is well known that the McVittie metric reduces to the 
Schwarzschild-de Sitter (Kottler) metric, which can be 
put into the static form pB)) . There is little doubt that 
in this case the central object described by this metric 
is a black hole. The Schwarzschild-de Sitter black hole 
has been the subject of much literature, mainly devoted 
to study the thermodynamical properties of dynamical 
horizons, which are interesting because of the simultane- 
ous presence of a Schwarzschild and a Rindler horizon. 
The inner horizon at f = to/2 has area 



This area does not depend on time because uih = 
for all McVittie solutions as a consequence of eq. PO)) . 



2to 



H 



(60) 



and is also time-independent. Therefore, the horizon 
area of the Schwarzschild-de Sitter black hole does not 
increase: the horizon does not expand, "resisting" the 
cosmic (accelerated) expansion. Cosmic censorship (in 
its version for non- asymptotically flat spaces 31|) is not 
violated in this case; the central singularity remains for- 
ever and is always surrounded by a horizon. 

It is not true that all physical systems in any expand- 
ing universe participate in the cosmic expansion; in a 
de Sitter universe this is true only for weakly coupled 
systems. This result agrees with the "all or nothing" be- 
haviour discovered by Price \2^. In the next section we 
argue that this phenomenon, however, is peculiar to the 
de Sitter background. 



BLACK HOLES IN ARBITRARY FLRW 
BACKGROUNDS: OLD AND NEW SOLUTIONS 

Although the Schwarzschild-de Sitter black hole does 
not participate in the cosmic expansion, in more general 
FLRW spacetimes black holes or other strongly gravitat- 
ing objects can expand under suitable conditions, as is 
shown below. 

The first step toward this discussion consists of iden- 
tifying exact inhomogeneous solutions of the Einstein 
equations that are suitable for describing strongly grav- 
itating objects embedded in a FLRW universe that is 
not a de Sitter space. In this case the McVittie met- 
ric (l39l) can not be interpreted as describing black holes 
[lol.]l3. TBI 55 1, however, it does describe some kind of 
singular central object. The computation of the area of 



the singular surface f — m/2 and of its curvature radius 
proceeds exactly as in the Schwarzschild-de Sitter case 
leading again to the time-independent area (j59p and ra- 
dius ([60]) . This applies also to the Reissner-Nordstrom 
generalization of McVittie's metric found by Gao and 
Zhang [l^. The issue is whether this central object is a 
realistic one. It is tempting to interpret the singularity at 
f = to/2 as an artificial one created by explicitly forbid- 
ding accretion of the cosmic fluid onto the central object, 
much like the axial singularity in cylindrically symmet- 
ric Bach-Weyl solutions [63|. In the latter, two massive 
particle are in static equilibrium at a finite distance and 
an axial singularity is interpreted as a strut holding them 
apart 6l|. One could think that the f ~ m/2 surface in 



the McVittie metric similarly acts as a wall to keep out 
the cosmic fluid. This interpretation seems to be corrob- 
orated by the following observation: the accretion rate 
for spherical accretion of a test fluid with energy density 
p and pressure P onto a Schwarzschild black hole of mass 
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/i is \6z 



r = TOo/2 is 



/i = AnDfi"^ {Poo + Poo) , 



(61) 



where poo and Poo are the energy density and pressure at 
spatial infinity, respectively, and is a constant corre- 
sponding to a first integral of motion [sl, [g^] ■ The accre- 
tion rate (I6ip vanishes for a fluid satisfying the quantum 
vacuum equation of state P — —p ~ Extrapolat- 
ing this result to a gravitating fluid, (spherical) accretion 
onto a Schwarzschild-de Sitter black hole is automati- 
cally prevented, the hole's mass does not change, and the 
condition (j40|) is satisfied naturally — it does no longer 
enforce the presence of a spherical wall to stop cosmic 
matter from falling onto the black hole. 

This interpretation of the 2-sphere singularity turns 
out to be at least partially erroneous. Below, some new 
exact solutions are presented in which the no-accretion 
condition is removed but the pressure is still singular on 
the surface r = m/2 unless the black hole is exactly co- 
moving. 



The Sultana-Dyer solution 

Recently Sultana and Dyer [l^ found an exact solu- 
tion of Petrov type D describing a black hole embed- 
ded in a spatially flat FLRW universe with scale factor 
a{t) cx t^/^ (in comoving time). The technique used to 
generate this solution consists of conformally transform- 
ing the Schwarzschild metric g^^'' — > 17^ g\^f^ with the goal 
of changing the Schwarzschild global timelike Killing field 

into a conformal Killing field for Sf^V 7^ 0, gener- 
ating the conformal Killing horizon (which differs from 
Hayward's future outer trapping horizons [63| and from 
dynamical horizons [g^l ) . The metric obtained by choos- 
ing f2 = a{t) = ao t^^^ — the scale factor of a dust-filled 
k = FLRW universe — is 



ds^ = -- 



2f 



(62) 

contrary to Ref. [iSII we use isotropic radius r and co- 
moving time t. This is formally the same as the McVit- 
tie solution ((39)) but with the important difference that 
the metric coefficient m is now a constant mo. To relate 
to our previous discussion, this implies that the Hawking 
quasi-local mass is mnit) = niQ a(t) and it increases with 
time. The condition is violated and an accretion flow 
of cosmic fluid onto the central object is present and is 
responsible for the increase in the gravitating mass. 

The source for this metric is a combination of two 
non-interacting fluids: Tat = T^^l^ + T^"\ where T^^' ^ 
pua Ub describes an ordinary (massive) dust and T^l^^ — 
Pn ka kb describes a null dust with density p„ and k'^kc = 
[3| ■ The surface area of the conformal Killing horizon 



Ait) 



dOdtp ^/gr. = 167rmQ a^(t) 



(63) 



and its physical (curvature coordinate) radius is simply 



2 Too 0-{t) J 



(64) 



which coincides with the familiar Schwarzschild radius 
r — 2m multiplied by the scale factor, as customary in 
FLRW cosmology. Thus, the physical radius of the hori- 
zon and the quasi-local mass are comoving with the cos- 
mic background. 

The Sultana-Dyer solution is non-singular at the sur- 
face f — m/2 and can be interpreted as describing a black 
hole embedded in a two-fluid universe. We see that re- 
moving the McVittie no-accretion condition can in- 
deed remove the singularity here, allowing the black hole 
to become comoving with the rest of the universe. There 
are, however, problems with the Sultana-Dyer solution: 
the cosmological fluid becomes tachyonic (negative en- 
ergy density) at late times near the horizon [l9j. More- 
over, it is desirable to have cosmological matter described 
by a single fluid composed of particles following timelike 
geodesies, and to drop the restriction to the special choice 
of the scale factor a oc t^/'^. 



McClure and Dyer [55| found another solution for a 
black hole-like object embedded in a radiation-dominated 
universe with a heat current, which satisfles the energy 
conditions everywhere and is perfectly comoving. How- 
ever, the energy density and pressure are singular at 
f = to/2, a similar solution in a dust-dominated uni- 
verse has singular energy density [5^. 

It is not clear at this point whether solutions exist 
describing black holes in arbitrary FLRW backgrounds, 
which are free of singularities at f = to/2 and satisfy the 
energy conditions everywhere. This is the subject of the 
next subsection. 



New exact solutions 

We look for solutions described by a generalized 
McVittie metric of the form ([39]) . but without impos- 
ing the no-accretion restriction (j40|) . and with arbitrary 
scale factor a(t). The line element is written in the form 



ds' = - 



where 



(t, f) 

^2 (t, f) 



dt^ + a^{t)A* (t, f) (rff2 -I- f^dn^) , 



A{t,f) = l + 



m{t) 
2f 



(65) 



S(t,f) = l-^. (66) 
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The only non-vanishing components of the mixed Ein- 
stein tensor are 



Gl 



3^2 
2m 



m 
fA 



m 

f^a'^A^B \ m 



(67) 
(68) 



GZ y^d 
9 — Lzo 



A2 



1 










2m 




fAB 



m 
fA 



It is convenient to introduce the quantity 

friH m B 



C: 



a 
a 



m 
fA 



rriH 



A 



(69) 



(70) 



which reduces to fnH/mH for Sultana-Dyer- type solu- 
tions with m ^constant. For any choice of the function 
771 (i) the quantity C reduces to 



Cy 



777 
777 



rriH 



(71) 



on the surface f = m/2. McVittie solutions have Ce = 
while Sultana-Dyer-type solutions have C = C^. — H 
everywhere. 

The Ricci curvature is 



R 



-g ^ab 



3^2 

"B2" 



2C + 



27)7(7 \ 

fAB J 



(72) 



and is singular on the surface f = 77i/2 if the pressure 
is singular there. We now specialize the discussion to 
different forms of matter. 



Perfect fluid 

If we assume that matter is described by a single per- 
fect fluid with stress energy tensor of the form 



Tab ={P + p) UaU(, + Pgab 



(73) 



and allow for a radial energy flow described by the fluid 
four- velocity u'^ = 0,0), it is easy to see that 

the only possible solution is the Schwarzschild-de Sitter 
black hole already considered. In fact, the normalization 
u'^Uc = — 1 yields 



A 



(74) 



and the Einstein equations give, using eqs. ([67|) - ((69|) 

niH = -GB^au (P + p) A\/ \ + a'^A^u'^ , (75) 



where A = J J dOdip ^fgr. = 4:^0^ A'^r'^ is the area of a 
spherical surface of isotropic radius f, 

^(^) = S'tG [(F + p) a^^V + p] , (76) 



= SttG [{P + p) a^A^u^ -f P] 



2C + 3C2 



2mG\ 
fAB J 



(77) 

= SnGP . (78) 



By adding eqs. ([77]) and ([78]) one obtains P = —p, i.e., 
only the de Sitter equation of state is allowed. Eq. ([75)1 
accordingly yields fnn = 0. With the exception of the 
non-accreting Schwarzschild-de Sitter black hole, a single 
perfect fluid can not source solutions representing spher- 
ically symmetric black holes embedded in a cosmological 
background. A mixture of two perfect fluids still consti- 
tutes a potential source, as exemplified by the Sultana- 
Dyer solution [l^. 



Imperfect fluid and no radial mass flow 

The following solutions describe perfectly comoving 
black holes. We assume now that cosmological matter 
is described by the imperfect fluid stress-energy tensor 



Tab = {P + p) UaUb + Pgab + QaUb + qb^a 



(79) 



where the purely spatial vector q'^ describes a radial en- 
ergy flow, 



,0,0,0 



(0,g,0,0) , (7=77^ = 



(80) 



and u'^Uc = —1. The (0,1) component of the Einstein 
equations Gab — SirGTab yields 

777 a AttG 

777 a 



■f^a^A^B^q 



(81) 



Since the Hawking mass is mn = m{t)a(t), it is 



rriH 



777 a 
777 a 



(82) 



and the area of a spherical surface of isotropic radius f is 
A = J J dOdip ^JgY. — 4:7Ta'^A*f'^, yielding the relation be- 
tween energy flow, area A, and accretion rate (Hawking 
mass added per unit time) 



mnit) = -GaB^Aq . 



(83) 



On a sphere of radius r >> m this can be written as (tak- 
ing into account the fact that radial inflow corresponds 
to g < 0) 



friH — GaA \q\ 



(84) 
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Hence, the quasi-local mass increases for inflow of matter. 
The (0,0) and (1,1) (or (2,2) or (3,3)) components 



J 



1 3^2 



SttG 
-1 



SttG S2 



a 


m ^ 


- 4 
a 


fAj 


d 


(- + 


di 


\ a 



rn 
fA 



of the Einstein equations yield the energy density and 
pressure, respectively. 



m 
fA 



31- + ^ 

a rA 



2m 
fAB 



(85) 



(86) 



from which it is clear that the energy density is always 
non-negative. The expansion scalar is 3C and, in terms 
of this quantity, eq. (I86p becomes the generalized Ray- 
chaudhuri equation 



and 



C 



In the limit 



3C2 
2 



> this reduces to the well known Ray- 

chaudhuri equation of FLRW cosmology H = 

AttGP, for which the Hamiltonian constraint = 
8ttGp/3 then yields 



H ^ -AttG (P + p) . 
Similarly, in the case m ^ 0, eq. ([55)1 yields 

fnC 

C^-^.G^iP^p)-^, 



(88) 



(89) 



which generalizes eq. 

It can be noted that, due to the factor B in the de- 
nominators, the energy density and pressure (|85p and 
(I86p appear to be singular on the surface f = m/2 where 
B vanishes. The situation is ameliorated by using the 
proper time r defined by dr = dt instead of the co- 
moving time t to absorb a factor B. This corresponds 
to using proper time instead of coordinate Schwarzschild 
time to offset an infinite redshift on the horizon of an 
ordinary Schwarzschild black hole, and yields 



A 
B 



m 
fA 



a 



rrir 
m 



as r 



so that 



SttGp = 3 



V a 



--A 



m/2 , (90) 



(91) 



8.GP = -2f^ + ^) -3f^ + ^ 
V a rAJ T \ a rA 



(92) 



The pressure, the energy density, and the Ricci scalar 
R — ~p + 3P appear to be finite on the surface f = m/2. 
This is true for any form of the function m(t) and con- 
trasts with the singularities in the solutions of Ref . |55| . 

Sultana-Dyer-type solutions with m = have a confor- 
mal Killing horizon describing a cosmological black hole, 
as in the Sultana-Dyer [l3| solution. By design, this black 
hole is perfectly comoving: also in this case the cosmo- 
logical expansion "wins" over the local strong field of the 
black hole. To keep friH > at f = m/2, where B — 0, 
one needs q ~* ~oo there. The Sultana-Dyer solution 
also suffers from a similar problem [l9| . This unphysical 
situation is due to the unrealistically simplified model of 
accretion. 



Imperfect fluid and radial mass flow 

We now consider an imperfect fluid with stress-energy 
tensor of the form (I79p and both radial mass flow and 
energy current described by 

u"" ^ (^^y^l + a^AH^, u, 0, 0^ , = (0, q, 0, 0) . 

(93) 

By using the components ((67)) -((69 l) of the Einstein tensor, 
the field equations become 



mn = -GaB^A\/l + a?A^v? [{P + p)u + q] 



(94) 
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-8ttG [{P + p) a^A^u" + p] 



(95) 



B I \ rAB 



= 8nG [{P + p) a^A*u^ + P + 2a^A^qu] , 



(96) 



(97) 



Adding the last two equations yields 
q = ^{P + p)-, 



(98) 



i.e., to an ingoing radial flow of mass there corresponds an 
outgoing radial heat current if P > —p. By substituting 
eq. (|98|) into eq. ((94l) . one obtains the accretion rate 



G 



niH = - — aSVl + aM'*u2 (P + p) Au , 



(99) 



where (P + p) Au can be seen as the flux of gravitating 
energy through the surface of area A. Since u < 0, the 
mass mn increases if P + p > 0, stays constant in a 
de Sitter background, and decreases if phantom energy 
with P < —p is accreted. This lends support to the 
conclusions of Ref. [s^j on the fate of a black hole in a 
phantom-dominated universe. 

Moreover, the energy density is given by 



A2 

8nGp= — 



3C^ 



G 



iC 



fAB 1 



'-A^i 



(100) 



For Sultana-Dyer-type solutions with m — mg ^constant 
the energy density reduces to 



8nGp 



P2 



2Ha^A^u^ 



1 



(101) 



and is positive-deflnite in a superaccelerating universe 
with H > 0, which is necessarily phantom-dominated 
[6^. Moreover, one can solve for the velocity of the fluid 
obtaining 



2a2y44 



1/2 



(102) 



The motion of the fluid becomes superluminal as f ^ 
mo/2, where B — > 0. In a realistic model the flow be- 
comes supersonic at a certain distance from this surface. 
This fact can only be taken into account in a more real- 
istic model of accretion, which will be studied elsewhere. 



DISCUSSION AND CONCLUSIONS 

It is well known that the effect of the cosmological 
expansion on weakly gravitating Newtonian systems of 
small size r (i.e., rHo << 1) is completely negligible for 
practical purposes, even though these systems do partic- 
ipate in the expansion. In larger systems such as voids, 
filaments, and large scale structures, the cosmic expan- 
sion plays a significant role BSlHlii. But the size 

of the local system is not the only relevant factor. The 
study by Price [22] is the first to focus on the strength 
with which the local system is bound. The "all or noth- 
ing" behaviour discovered constitutes an important step 
in the understanding of the process. However, Price's 
discussion is limited to a de Sitter background which is 
too special to draw general conclusions. Moreover, the 
classical non-relativistic "atom" considered can not de- 
scribe arbitrarily strong binding of the local system be- 
cause, when the energy density and stresses involved be- 
come very large, they distort spacetime causing the met- 
ric to substantially deviate from a cosmological one. It 
is preferable to study exact solutions describing a local 
object with a strong local gravitational field — e.g., a 
black hole — embedded in an expanding universe. 

Independent motivation for our study arises from the 
recent realization that if the current acceleration of the 
universe is dominated by phantom dark energy with 
P < — p, then the universe r nay be running into a Big Rip 
at a finite time in the future [30| . Recent literature has fo- 
cused on the way this catastrophic accelerated expansion 
of the universe comes to dominate the local dynamics of 
bound systems (clusters, galaxies, stars, etc.) and tears 
these systems apart. In this context, it is interesting to 
pose the question of whether the Big Rip can destroy 
a black hole horizon and expose the central singularity, 
thus violating cosmic censorship [sij. A partial answer 
comes from Ref. 32| : these authors analyse spherical ac- 
cretion of a phantom test fluid onto a Schwarzschild black 
hole and, extrapolating the results to a gravitating fluid, 
reach the conclusion that the horizon disappears before 
the Big Rip together with the central singularity, without 
violating cosmic censorship. This result is quite plausi- 
ble, however it is desirable to have a study of the accre- 
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tion process for a gravitating fluid, which brings us again 
to the realm of exact solutions. 

The fully relativistic systems considered in the present 
paper provide at least a partial answer to the ques- 
tions above. The McVittie solution (|39l) and (|40l) [ij 
is accretion-free, describes a general FLRW background 
universe, and does not expand. However, it has a mild 
singula rit y on the surface f = m/2 (the putative horizon) 
[l0|,llj,|l5| and therefore it does not describe a black hole. 
There is an important exception, the Schwarzschild-de 
Sitter black hole (a special case of the McVittie metric) 
which does not suffer from this singularity and describes 
a true black hole horizon which resists the cosmic expan- 
sion. This feature is consistent with Price's study of the 
classical atom in a de Sitter background [13] and with 
our phase space picture of Sec. 2.3. 

The overall picture emerging is that the considera- 
tion of a de Sitter background is rather misleading, even 
though it considerably simplifies the calculations: the 
"all or nothing" behaviour is not generic of FLRW space 
but is limited to a de Sitter background. The next sys- 
tem considered, the Nolan interior solution [s^ , does not 
suffer from the singularity problem of McVittie's metric 
because the surface f = m/2 is covered by the matter 
composing the star. This solution, which can be thought 
of as providing a source for McVittie's metric to which it 
is matched, does not accrete either and can be embedded 
in a general FLRW background: it is comoving. Remov- 
ing the limiting assumption of a de Sitter background 
allows the strongly gravitating central object to expand. 

Another exact solution describing a black hole embed- 
ded in a non-de Sitter universe is the one of Sultana and 
Dyer [l^ in which, contrary to the McVittie metric, ac- 
cretion onto the black hole does occur and the confor- 
mal Killing horizon (the black hole horizon) is, by de- 
sign, comoving. The peculiarity of the de Sitter cosmos 
is thus further put into evidence. The Sultana-Dyer so- 
lution, however, suffers from the following limitations: 
i) it is restricted to the special form of the scale factor 
a{t) = aQt^/^ in comoving time; ii) the matter source is 
not a simple fluid but a mixture of two non-interacting 
fluids, one of which is a null dust; and iii) the energy 
density becomes negative at late times near the horizon. 
It is well known that, in general, matching black hole 
and cosmological metric produces stress-energy tensors 
that violate the energy conditions in some regions of the 
spacetime manifold |55j . 

To overcome these difficulties we have studied new al- 
ternative solutions of the Einstein equations which gen- 
eralize the McVittie metric by allowing radial accretion 
onto the central object. We have presented new solu- 
tions for which the surface f = m/2 is non-singular and 
is perfectly comoving. These new solutions do not always 
suffer from the limitations i)-iii) above — in particular, 
the energy density is everywhere positive at all times for 
one of the solutions, but the accretion flow generally be- 



comes superluminal, an artifact of the simplified "rigid" 
model of accretion used, in which matter everywhere in 
the universe moves toward the central black hole (albeit 
its radial speed becomes zero far away from it). More- 
over, these solutions support the result of Ref. [s^l that 
a black hole embedded in a phantom-dominated universe 
disappears, respecting cosmic censorship. 

It appears, therefore, that the strong external cosmo- 
logical gravitational field can distort a black hole hori- 
zon, or anyway stretch an object dominated by a strong 
local field, de Sitter-like expansion does not, but this 
should be seen as an exception to the rule due to the 
scale-invariant nature of the exponential scale factor (in 
fact, the de Sitter metric can be put into static form). 
A posteriori^ the fact that a black hole horizon expands 
with the cosmic substratum should not be regarded as 
surprising: many exact solutions are known in which 
a black hole horizon is distorted by its surroundings, 
due to an external gravitational or acceleration field 



11, 37, 38 



39 



41, 42, 43, kj, an electric M 



^iii, iiiii, 1^,^. i^jyi. i^^i. 1^^. ntj. . i^TK Liiaj, t^ii ^^^^^^^^ or 

mag netic [42[13l4al5^ field, or combinations of them 
5ll | . A substantial amount of literature has been devoted 
to horizon deformation ( [s^ and references therein). 

Our results do not constitute the last word on the is- 
sue of cosmological expansion and local systems. Future 
endeavours include the search for more general and more 
realistic exact solutions of the Einstein equations describ- 
ing accreting black holes in arbitrary FLRW or Bianchi 
backgrounds. 
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